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Abstract 

i-C ■ We consider for a monatomic liquid the density and current autocorrelation functions from the 

point of view of the Vibration- Transit (V-T) theory of liquid dynamics. We also consider their 
Fourier transforms, one of which is measured by X-ray and neutron scattering. In this description, 
the motion of atoms in the liquid is divided into vibrations in a single characteristic potential 
valley, called a random valley, and nearly-instantaneous transitions called transits between valleys. 
The theory proposes a Hamiltonian for the vibrational motion, to be corrected to take transits 
into account; this Hamiltonian is used to calculate the autocorrelation functions, giving what we 
call their vibrational contributions. We discuss the multimode expansions of the autocorrelation 
functions, which provide a physically helpful picture of the decay of fluctuations in terms of n-mode 
scattering processes; we also note that the calculation and Fourier transform of the multimode series 
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are numerically problematic, as successive terms require larger sums and carry higher powers of the 
temperature, which is a concern for the liquid whose temperature is bounded from below by melt. 
We suggest that these problems are avoided by directly computing the autocorrelation functions, 
for which we provide straightforward formulas, and Fourier transforming them numerically. 
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I. INTRODUCTION 



V-T theory is being developed specifically to describe atomic motion in real monatomic 
liquids. This theory adopts a point of view different from, and in a certain way comple- 
mentary to, the more traditional formulations. The theory most successful for equilibrium 
thermodynamic properties is based on the interatomic pair potential and the pair distri- 



bution function 



2j. This theory can be made exact by including n-body potentials for 
n > 3. On the other hand, the most widely studied formulation of nonequilibrium proper- 
ties, as expressed through time correlation functions, is generalized hydrodynamics pi, U, [ij] . 
The memory function analysis is formally exact, but has not yet proved amenable to first- 
principles evaluation. V-T theory is based on an approximate but tractable zeroth order 
Hamiltonian, which in principle can be systematically improved to include corrections be- 
yond zeroth order. It complements traditional theory in that both equilibrium and nonequi- 
librium properties can be calculated from the same zeroth order theory, without adjustable 
parameters 

V-T theory is based on a symmetry classification of the potential energy valleys among 
which the liquid atoms move jf], B, [^j] • Specifically, the amorphous (noncrystalline) valleys 
are classified as (a) symmetric valleys, with e.g. microcrystalline structures, which have 
significant variations in their potential energy and vibrational frequencies but are relatively 
few in number, and (b) random valleys, with approximately random close-packed structures, 
which all have the same potential energy and vibrational frequencies in the thermodynamic 
limit, and which are of overwhelming numerical superiority. As a zeroth order model for 
the potential energy surface, we neglect the symmetric valleys, and construct an "extended 
random valley" as the harmonic extension to infinity of one (any) random valley. The 
partition function is then the extended random valley partition function, times the number 
of such valleys. This is a universal number, calibrated from the constant density melting 
entropy of the elements [s| . The zeroth order theory gives a good account of the experimental 
thermodynamic properties of the elemental liquids at melt, indicating that the corrections 
are small at melt Q]. 

In a real monatomic liquid the motion consists of two parts: vibrational motion in a 
random valley, and transit motion, wherein the system crosses valley- valley intersections. In 
applying V-T theory to a physical process, we evaluate the vibrational contribution in zeroth 
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order, exactly and without adjustable parameters, and then model the effect of transits. 
Thus in the liquid partition function the major contribution is vibrational, while the role 
of transits is to allow the system to access the multitude of random valleys, and hence 
to possess the extra entropy associated with melting. In the dynamic response function we 
have found that both the location and area of the Brillouin peak are given by the vibrational 
contribution, while transits are responsible only for a partial broadening of this peak. Hence 
an a priori calculation of the vibrational contribution to the Brillouin peak dispersion curve 
for pseudopotential sodium agrees within experimental error with the measured results of 
inelastic x-ray scattering for real liquid sodium [h]] . In addition, the vibrational contribution 
can be modified to account for transits, and this model provides an excellent fit to the entire 



Brillouin peak 



Hi 



The purpose of the present paper is to study the vibrational contributions to density and 
current correlations in a monatomic liquid. The density autocorrelation function F(q, t) 
and its Fourier transform S(q, u) are the basis for inelastic neutron scattering theory. This 
theory was originally developed for crystals, and expansions of S(q , u) in powers of nuclear 
d^acements were a„a.yzed in detail by Maxadndin and Fein Q, Ambegaokar, Conway, 
and Baym [13], Cowley [14] . Lovesey [15|, and Glyde [16|. A comparison of different ex- 
pansions may be found in Glyde 16(. (Physically motivated approximations for the glassy 
state are discussed by Egelstaff [17].) Liquid theory has important similarities to, and dif- 
ferences from, the crystal theory. First, in liquid and crystal alike, the multimode expansion 
converges more slowly as the momentum transfer q increases, and as the temperature T in- 
creases. While the range of q is approximately the same for inelastic scattering from a given 
material in crystal and liquid form, the temperature is quite the opposite, being below T m for 
the crystal and above T m for the liquid. Second, vibrational dynamics is more complicated 
for a random valley, because every site is inequivalent, as opposed to inequivalence of only 
the crystal sites in one unit cell. As a result, the dynamical matrix is 3N x 3jV for a random 
valley, but is block diagonal in N 3 x 3 matrices for a primitive lattice crystal. We will show 
how to overcome convergence problems resulting from the high temperatures of the liquid, 
and how to address computability problems resulting from the large dynamical matrix. On 
the other hand, we will not address anharmonic corrections here, as it is first necessary to 
formulate and test the harmonic theory. Further, since the atomic motion is classical in 
most monatomic liquids, we work with classical statistical mechanics, which allows us to 
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freely commute atomic positions and momenta. 

The longitudinal and transverse current autocorrelation functions CV t(q, t) are discussed 

n n n 

by March [2|, Hansen and McDonald 4J, and Balucani and Zoppi [5j. The relation of these 
functions via Green-Kubo theory to the longitudinal and transverse viscosities is described, 
as well as their treatment within the framework of generalized hydrodynamics. We no- 
tice that Ct(g, lo) = (uj 2 /q 2 )S(q,uj), while CV(<?, co>) is an independent function 0,0]. For 
the liquid alkali metals, dispersion curves associated with C^g, to) have been attributed to 
sound-like collective modes of the atomic motion (|5|], Fig. 6.6). These authors also argue for 
the existence of distinct fast and slow decay mechanisms in Ci{q,uj) for monatomic liquids 
( 5I, Sec. 6.2). In V-T theory, these would result from vibrations and transits respectively. 
Dispersion curves associated with Cx^q^oj) have also been studied [l8[]. For the transverse 
current correlation function in liquid rubidium, it has been shown that the usual phenomeno- 
logical relaxation time approximations for the memory functions are not capable of fitting 



the MD data 



19]. 



Currently there is no theory, short of computer simulations, that will account for these 
dispersion curves without adjustable parameters. Since the vibrational contribution alone 
does this for the Brillouin peak dispersion curve [10], it is possible that it will do so as well 
for the dispersion curves associated with C^q, oj) and Criq, oj). This possibility motivates 
the present work. In addition, the computational problems mentioned above for S(q,u) 
need to be addressed also for Cl,t(<1, oj); these are the poor convergence of the multimode 
expansions at high temperatures, and the difficulty of evaluating nested 3iV-fold sums for 
moderately large N. 

In Sec. [Til we briefly summarize the vibrational normal mode theory for an extended 
random valley, and obtain the pair correlation functions involving atomic positions and 
velocities. In Sec. IIHI we write the formulas for vibrational contributions to the dynamic 
response functions, F(q, t) and S(q,u); we compare these with the formulas derived by 
Carpenter and Pelizzari 20J]; and we show how to circumvent practical limitations posed 
by convergence and computability problems. The vibrational contributions to Ci,T{q,t) are 
derived in closed form in Sec. IIVI and are shown to consist of decoupled velocity-position 
terms and velocity-position interference terms. The one- and two-mode contributions to 
Cl,t(<1,u) are derived, and accurate but practical evaluation procedures are described for 
all the current correlation functions. A brief summary is presented in Sec. [V] 
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II. RANDOM VALLEY VIBRATIONAL MODES 



This topic has been discussed previously [211 ]. so a brief outline will suffice here and will 



establish our notation. Our mechanical system consists of N atoms in a cubical box, with 
the atomic motion governed by periodic boundary conditions at the box surface. The atomic 
positions and momenta at time t are respectively rx{t) and Pi<(t) for K = 1, . . . , N. For 
motion in a single random valley we follow the lattice dynamics tradition of Born and Huang 



22| and write 

r K {t) = R K + u K {t), (1) 

where Rk is the equilibrium position and WR-(i) is the displacement. The system potential 
is $({r^}) and the potential at the equilibrium structure is $o = The extended 

random valley is by definition harmonic, and has the Hamiltonian 

p 2 1 

H mb = % + J2^ + 2^^ ^Ki, Lj UKiU Lj , (2) 
Ki KL ij 

where i,j represent Cartesian directions and M is the atomic mass. The potential energy 
coefficients &Ki,Lj form a real symmetric 3N x 3iV matrix, which is the dynamical matrix 



22l | for the random valley. The vibrational normal modes have coordinates and momenta 
respectively q\(t) and p\(t), for A = 1, . . . ,3N. The transformation from displacements is 
defined by 

qx = WKi,\ UKi, 
Ki 

PX = ^W K i,XPKi, (3) 
Ki 

where the eigenvector components Wxi,x are real and diagonalize the dynamical matrix. It 
follows that 



2M 2 



Ti-vib = $0 + 

A 

where uj\ are the mode frequencies, from which it also follows that 



(4) 



qx{t) = qx{0)cosuJxt + ^-^-smuJxt 

Px{t) = Pa(0) cosuxt - Mu\q\(0) sincc^t (5) 
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The eigenvector components also satisfy orthonormality and completeness relations 21], 
respectively given by 

y^w KLX w K i,\> = 6xx>, (6) 

Ki 

y^ J WKi,\WLj,\ = S K LSjj. (7) 
A 

The basic ingredients of vibrational time correlation functions are the mode-mode func- 
tions such as {qx(t)qx'(0))- These can be calculated from Eqs. (j4j) and ((Sj), in classical or 
quantum mechanics, where the classical limit of the quantum result is the leading term in 
the expansion for kT/hwx > 1. For — oo < t < oo, the classical expressions are 

kT 

(q\(t)q\>(0)) = Tf^ S ™ cos ^, (8) 

A 

(P\(t) PA'(O)) = MkT8xx> cos uj x t, (9) 
kT 

{px(t)qx>(0)) = 5 xy smu x t. (10) 

^A 

Correlation functions involving the atomic positions and momenta follow from these and 
the normal mode transformation, Eqs. ([3]). It will be convenient to use velocity Uk instead 
of px = Muk, and to eliminate the Cartesian components in favor of vector dot products. 
Then with a and 6 arbitrary vectors, 

kT \ - 

(a ■ u K (t) b ■ u L (0)) = — 22 a ' w k\ b ■ w lx cosuj x t, (11) 

A 

/ / \ i / w kT , sino;^ / s 

(a ■ u K [t) b ■ u L (0)} = — — > a ■ w KX b ■ w LX , (12) 

M LOx 

i / \ i i w kT ^-^ cos ujxt , \ 

{a ■ u K (t) b ■ Ul(0)) = — > a ■ w KX b ■ w LX — . (13) 

M u x 

The fluctuation results are obtained by setting t = 0. Notice Eq. (JT2]) . the velocity- 
displacement correlation function, is odd in t, and the fluctuation vanishes. Eq. ( Tl2|) satisfies 

(a ■ u K {t) b ■ u L (0)) = -(a- u K (t) b ■ u L (Q)) , (14) 

which is valid for any motion as a consequence of time translation invariance. Eq. ( 1121) also 
satisfies 

(a ■ u K (0) b ■ u L (t)) = (a ■ u K (t) b ■ u L {0)) , (15) 
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but this relation follows from the specific form of Eq. (|T2|) and is not true for general 
motions. The above equations are valid at all temperatures for which the atomic motion 
is classical. Hence these equations are appropriate for comparison with ordinary (classical) 
MD calculations at all T > 0. 



III. DYNAMIC STRUCTURE FACTOR 



Because the system motion is constrained by periodic boundary conditions, Fourier com- 
ponents are defined only at a discrete set of q. To reflect the isotropic nature of a real 
liquid, a macroscopic average can include an average over the star of q for each \q\, this 
average being denoted (•••)„»• The density autocorrelation function, also called the van 
Hove function, is then written 

F(g,t) = i^X) e ^ ( * ) E e< * rjl(0) )) • ( 16 ) 

Here each sum is a q-component of the system density, and the inner average is over the 
vibrational motion, as in Eqs. (IHl fToT) . Notice the q* average makes F(q,t) a function only 
of the magnitude q = \q\. We have previously used notation such as F vib (q,t) and (• • - ) vib , 
but here the vibrational subscript will be omitted. Then inserting Eq. ([TJ for r^(t), 

F(q,t) = ^(^e^ R « L (e^< u *®^ »)\ , (17) 



KL 



where Rrl = Rr — Rl- The vibrational average is evaluated by Bloch's identity (e.g. [15]), 
which is derived via the normal mode transformation, to write 



F(q,t) = — /^ e -^-^ e -f ([9-("K(i)-™z,(0))] 2 ) \ 

^ \KL la* 



(18) 



Finally we separate the fluctuation terms to find 

\KL I q * 

where e~ Wli<yC ^ is the Debye- Waller factor for atom K, which can be expressed with the aid 
of Eq. (PSD as 

^( 9 H§£^f< (20) 



The function in the exponential may also be evaluated using Eq. ({TBI . 

The time dependence of F(q, t) is contained in the time correlation functions 
(q ■ u K (t) q ■ u L (0)) in Eq. (|TU|) . The multimode expansion is obtained by expanding in 
powers of these functions, according to 

e ( q . UK (t) q .uU0)) = l + (q. UK (t) q ■ u L (0)) + ~{q- u K (t) q ■ u L (0)) 2 + ■■■ . (21) 

This is a useful procedure, because each term can be Fourier transformed analytically, 
to give an explicit picture of the scattering process in each order. For a random valley, 
(q ■ u K {t) q • u L (0)) — > as t — > oo, so that F(q, oo) is positive and is given by 

2\ 



F(q, oo) 



Then the series for S(q,u) is 



N 



K 



-iq-R K p -W K {q) 



(22) 



S(q, to) = F(q, oo)5(u) + S^{q, to) + S®(q, u) + ■ ■ ■ . 



(23) 



The leading term expresses elastic scattering. S^'(q,uj) arises from the linear term in 
Eq. (l2~Tj) . and with Eq. ([13"]) the expression can be organized into 



S {1 \q,u;) 



M N 



(24) 



where 



e -w K { q ) q . Wrx _ 



(25) 



(Note that this f\(q) differs from our definition in 10, although S^(q,uj) is the same.) 
Hence (q, u) is the sum over all vibrational modes of single-mode scattering at momentum 
transfer hq, in which energy hio\ is either created in mode A (the S(u + u\) term), or is 
annihilated in mode A (the S (u> — u\) term). 
In a similar way, S^(q,uj) can be written 

' kT\ 2 1 

^£(I/aa'(?)|%<%±^±o; a 0>s, (26) 



S (2 \q,co) 



M 2N 



where 



fxx'{q) 



-E 

din j—^ 



e - iq .R K e -w K{q) q . WKxq . Wkx , 



(27) 



K 
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and (8(uj ± u\ ± ojy)) s is the symmetric average of four (^-functions, 

(5(UJ ±UJ X ±UJ X ')) S = ~ 5(U + UJ X + UJy) + 5(U - UJ X - UJy) 

+5(cu + uj x - ujy) + 5(tu - lu x + ujy) . (28) 

is therefore the sum over all pairs of vibrational modes of the cross section for scattering 
at momentum transfer hq, in which two excitations are created (the 5(u + uo x + u x >) term) 
or annihilated (the 5(u — uj x — ujy) term), or in which one excitation is created and another 
is annihilated (the 5{uj ± (u x — uj x >)) terms). 

In the multimode expansion of S(q,to) we see from Eqs. ff24|) and fl26|) that S^ n '(q,u) is 
proportional to T n . Therefore any finite termination of the series diverges with increasing 
T even though S(q,cu) does not. We also see from the above equations that S^ n \q,u) 
requires evaluation of n + 1 nested sums, each over atoms or over normal modes, so that 
each sum contains of order of N terms. To characterize the magnitude of the computation, 
we consider a system of iV = 1000 atoms, for which the data {-R^; uj x , w KX } are available. 
Then single, double, and triple sums are readily computed, while a quadruple sum challenges 
current computer facilities. On the other hand, problems related to both convergence and 
computability of the multimode expansion are eliminated by writing 

F(q, t) = [F(q, t) - F(q, oo)] + F(q, oo) (29) 

where F(q, oo) is given in Eq. ff22l) ; the first term can be Fourier transformed numerically, 
and the transform of the second term is just F(q,oo)6(u) (see Eq. (|23|) ). Since the time 
correlation functions in Eq. ( [19]) are proportional to T, F(q,t) is well behaved as a function 
of T. Further, F(q,t) can be computed exactly as a triple sum, as can any term in the 
multimode expansion of F(q,t). 

The quantum expressions for S^Hq.cu) and S^ 2 \q,cu) for a polyatomic amorphous solid 



were derived by Carpenter and Pelizzari [20j. When their results are applied to a monatomic 
system and the classical limit is taken, our Eqs. (I24H28P are obtained with the following minor 
differences. Our eigenvectors are real (not complex) and are normalized to 1 (not N); our 
A' = A term in Eq. ( |26|) contains the two inelastic terms of Carpenter and Pelizzari, plus 
two additional elastic terms in 5(uj)] and by writing the Y2kl m ^1- CG3) a s | YIk' " \ 2 
in Eqs. ( I241T271) . our multimode terms have one lower order of nested sums than those of 
Carpenter and Pelizzari. 



IV. CURRENT AUTOCORRELATION FUNCTIONS 

Fourier components of the particle current are denoted j(q,t), where 

j( 9 ,t) = 5> J r(t)e- , '"-*e>. (30) 

K 

Let us define for each q the orthogonal set of unit vectors q, Si and s 2 . Then the longitudinal 
and transverse currents are 

j L (q,t) = qq j(q,t), 

j T (q,t) = • j(q,t) + s 2 s 2 - j(q,t). (31) 
The corresponding current autocorrelation functions are defined by [s| 

C L (q,t) = ±((j L ( q , t ).f L (-q,0))) q *, 

C T (q,t) = 2^«Jr(flf,t)-i5.(-9,0)» tf .. (32) 

In these defining equations the inner averages are over the complete atomic motion for a 
macroscopic system. In what follows we shall limit ourselves to the vibrational motion in a 
finite system with periodic boundary conditions. 

Let a represent any one of q, 8±, s 2 , and define the autocorrelation function 

D(a, q, t) = i ( ( « • r K (t) e^"*® ^ a ■ r L (0) e^) \ \ , (33) 



A 



q 



where the inner average is over vibrations in an extended random valley. Then both Cl and 
Ct are represented by D, since 

C L (q,t) = D(q,q,t) 

C T (q,t) = l -[D{s liq ,t) + D(s 2 ,q,t)]= D(a 1>g ,t). (34) 

The last equality reflects the isotropic symmetry of a real liquid. 
With Eq. (CED for r K {t), Eq. (EHJ becomes 

D(a jq ,t) = ^(Y,e- iqRKL ((a-u K (t)a-u L (0)) e -^ K (t)-u L (o))A (35) 



The average is quadratic in velocities, and contains terms of all positive powers in displace- 
ments. The odd powers are imaginary, and harmonic averages involving them vanish. Hence 
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the nonzero averages are quadratic in velocities and of all even powers in displacements. To 
evaluate the thermal average in Eq. fl35l) . we first use the normal mode expansion from 
Eq. ([3]), with the result 



(a ■ u K {t) a ■ u L (0) e -*9-[«*(t)-«z.(o)]) 
= y^Q" ™K\ a ■ w LX > ( q\(t)q\'(0) exp 



A, A' 



(36) 



where the vectors Wkx have components WKi,x- Eliminating q x (t) and q x (t) using Eqs. (J5]) 
and replacing q\(0) and p x (0) by q x and p x respectively yields 

a • ^^Afl ■ w LX > ( ( C ° S A ^ 2 xt PxPy - UX ^ Uxt q\P\' ) exp ( ^ A n% + B tf>v ) ) ( 37 ) 



A, A' 



where the coefficients and can easily be calculated, but their exact values are un- 
needed. The first term in brackets is proportional to 



P\P\> 



exp A ^ + B ^t^j 



d d 



_d d_ 

dB x dBy 



exp 



dB x dB x 

2 



- ^exp A »% + B fPv J J 



(PxPx>) + ( \J2A u q u + B : 



vPv PX 



^2A u ,q v ,+B, 



v'Pv' Px> 



x exp 



\(^{Y, a ^ + b ^ 



[(j>\P\>) - (q ■ [u K (t) - u L (0)} px) (q ■ [u K (t) - u L (0)} px>) ] 



x exp 



--((q-[u K (t)-u L (0)]) 2 ) 



(38) 



In passing to the second line we used Bloch's identity, 



exp (^£2 4*9** + j J = ex P A M + j 



(39) 



(Lovesey's proof 15|] is for the quantum mechanical single-mode case; the generalization to 
multiple modes is trivial. The result can also be proved very simply in the classical case by 
completing the squares in the integrals.) In the third line we relabelled dummy indices, and 
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in the final line we returned the quantity in the exponential to its original form. A similar 
calculation shows that the other term can be expressed 

d d 



exp A v% + B ^Pv J 



dA x dBy 



exp 



[(q\p\>) - (q ■ [u K (t) - u L (0)\ q x ) (q ■ [u K (t) - u L (0)] p x >) ] 
-\{{q-[u K {t)-u L mf) 



x exp 

Putting this back together and rearranging, we get 
(a ■ u K (t) a ■ u L (0) e -w-[«*(*)-«*(°)]) 



(40) 



22 « • w kx a ■ w LX 

A, A' 



COS U\t 



COS U\t 

M 2 



-P\ - 



M 2 
ui\ sin uj\t 



u>x sin uj\t 

w 



M 



q\)q- [u K (t) - u L (0)] > (q ■ \u K {t) - u L {0)] p y ) 



x exp 



--((q-[u K (t)-u L (0)]f) 



(41) 



Next, we change q\ and p\ back to q\(Q) and p\(0) and use Eq. (jSJ) to simplify the time- 
dependent factors, reducing the expression to 



^ a ■ w KX a ■ w LX , 

A, A' 



(qx(t)qx>(0)) ~ (q ■ [u K (t) - u L (0)} q x (t)) 



x (q ■ [u K (t) - u L (0)) q x ,(0)) 



exp 



--((q-[u K (t)-u L (0)]) 2 ) 



(42) 



When this is summed over A and A' using Eq. ([3D, the expression becomes 



[(a ■ u K {t) a ■ u L (0)) - (q ■ [u K {t) - u L (0)} a ■ u K {t)) (q ■ [u K (t) - u L (0)) a ■ u L (0))] 



x exp 



1 



{q-[u K (t)-u L (0)}) 2 ) 



(43) 



The second term, involving products of displacements and velocities, can be simplified some- 
what. Terms such as (q ■ Ux(t) a ■ Ux(t)) are invariant under time translation, so they can 
be evaluated at t — 0. In this case, the expression is a linear combination of terms of the 
form (q x p x >), all of which vanish. Thus the two equal-time terms drop out, so 

(a ■ u K {t) a ■ u L (0) e -i^K(t)- UL (p)]\ 



[{a ■ u K (t) a ■ Ul(0)) + (q ■ u L (0) a ■ u K (t))(q ■ u K (t) a ■ u L (0))] 



x exp 



i<(g-[^(t)- WL (0)]) 2 > 



(44) 
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The exponential is the same term that appears in Eq. (fl8|) . so it may be transformed the 
same way. Putting all this together, Eq. fl35l) becomes 



D(a,q,t) = —I ^ e -iqRKL e ~W K (q) e -W L (q) e (q-u K (t) q-u L (0)) 



N 4 

\ KL 

[a-u K (t) a ■ ii L (0)) + {a ■ ii K (t) q ■ u L (0)) (a ■ ii L (0) q ■ u K (t)) \ , (45) 

/ q* 

and all of the correlation functions may be evaluated with the aid of Eqs. (11111131) . Compar- 
ison with Eq. ( Tl9i) verifies the well known relation for the longitudinal correlation function 

D(q,q,t) = -\F(q,t). (46) 

Just as with F(q,t), expansion of the time- dependent exponential in Eq. fj4"5l) produces 
a series which can be Fourier transformed analytically, term by term. We express this 
trasformed series as 

D(a, q, u) = D (1) (a, q, w) + D (2 \a, q, v) + - • • , (47) 

where the superscript reveals the number of normal modes coupled in the correlation. D^ 1 ' 
is the pure velocity term, given by 

£>«(&, q, ^) = ~£ q)\ 2 ) q , \ [8{u> + wa) + 8{u - w A )] , (48) 

A 

where 

fc A (a, <?) = ^ e - J 9 «x e -H/K( 9 ) & . (49) 
if 

The similarity with one mode scattering as described by Eqs. (l24"j) and (|25p is apparent. The 
disappearance here of the tu^ 1 factor in Eq. (I25I) arises from velocities replacing displace- 
ments. 

The two-mode contribution takes the form 

17 ) N^TW^'^^ Kx'&Q^q* {S(uJ±UJ X ±UJy)} s 



17 ) mJ2 ; (^'(a, <?) ^a'a(o> 9)),. ± ± , 

' AA' ^ A ^ A ' 



(50) 



where 

/i AV (a,q) = ^e-^'e-^^ a-^ A (51) 
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and where (8{u ± uo\ ± LUy)) A is the antisymmetric combination of four ^-functions, 

(5(lU ±LUx±LU X ')) A = ^ 6(co + LU X + LU X ') + 5(lj - LU X - U\>) 

-5(lu + lu x - lu x >) - 5(lu - LU X + LUx>) ■ (52) 

The first line of Eq. (!50|) is the velocity-decoupled term, and contains (\^2 K ■ ■ ■ \ 2 ) gt , similar 
to S( 2 '(q,u) in Eqs. fl26l) and fl27|) . The second line of Eq. fl50|) is the velocity-interference 
term, and does not contain an absolute square, because of the interchange of A and A' in h xx ,. 
It is still possible to show that the velocity-interference term is overall real, even without the 
q* average. A more prominent difference from the velocity-decoupled term is the change in 
sign of two of the (^-functions in the velocity-interference term. This will give the two terms 
a strongly different ^-dependence, since 5(lu ± (lux + ujy)) * s spread broadly over lu, while 
S(lu ± (lux — lu\')) is concentrated at smaller frequencies. 

V. SUMMARY 

In V-T theory of monatomic liquid dynamics, the vibrational contribution to a statisti- 
cal average can be evaluated exactly from the vibrational eigenvalues and eigenvectors of a 
single random valley. The transit contribution consists of small but not negligible correc- 
tions to the vibrational contribution, and has so far been accounted for by parametrized 
models. In the thermodynamic properties of the elemental liquids, the transit correction 



accounts for the constant density entropy o 
mal energy from valley-valley intersections 



meltine; and for the contribution to ther- 



231 ] . In nonequilibrium properties, as expressed 



by time correlation functions, the transits cause decay of correlations in addition to that 



24j . Continuing in this line of the- 



already present in the vibrational contribution 
oretical development, we have two goals in the present paper: (a) to derive the equations 
for the vibrational contribution to the current autocorrelation functions, and (b) to address 
difficulties present in the liquid state in computing these functions and the related density 
autocorrelation functions. 

The current autocorrelation functions are denoted D(a,q,t), where a = q for the lon- 
gitudinal current, and a is transverse to q for the transverse current. A closed form ex- 
pression for D(a,q,t) is derived in Eq. (T4"5]) . This expression contains two classes of terms: 
(a) velocity-decoupled terms, containing one velocity-velocity correlation function, and (b) 

14 



velocity-interference terms, containing a product of two velocity-displacement correlation 
functions. The one- and two-mode contributions to the Fourier transform D(a,q,u) are 
written in Eqs. (I48H52I) : D^'(a, q, u) expresses decay of current fluctuations through pro- 
cesses involving one vibrational mode, and D^ 2 \a,q,u) expresses the same decay through 
processes involving two vibrational modes (including one mode twice). 

The density autocorrelation function F(q,t), and its Fourier transform S(q, u>), are the 
basis of inelastic neutron scattering theory for condensed matter. With notable exceptions, 



the point of view in crystal theory 
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Q, and in amorphous solid theory 2(|, is 



that the multimode expansions of these functions converge well. The exceptions include 
ihe quantum crystals 16j|, and soft phonons associated with the hcp-bcc transition in Ti 



2a 



261 ] and Zr [23, |28(. However, this viewpoint is precarious for the liquid, because 



the characteristic liquid state persists to temperatures of several times T m |29j |. while the 
multimode expansions in question diverge in T at any finite order. The same convergence 
issue is present in the current correlations, and the following conclusions apply to density 
and current correlation functions alike. 

(a) The multimode series for the vibrational contribution to the Fourier transforms S(q, u) 
and D(a, q,oo) are very informative, as they resolve the decay of fluctuations into one-mode 
processes, two-mode processes, and so on. In classical statistical mechanics, the one-mode 
processes go as T and are computable as a double iV-fold sum; the two-mode processes go 
as T 2 and are computable as a triple iV-fold sum, and so on. In general we expect the series 
to converge at about the same rate for liquid and crystal phases at T m . Then, as T increases 
above T m , higher order terms which are more difficult to compute become more important. 

(b) These convergence and computability problems are circumvented by computing the 
vibrational contributions to the time correlation functions, F(q,t) and D(a,q,t), and by 
numerically transforming these. (We suggest subtracting off F(q, oo) and transforming that 
term analytically, which is not necessary for D(a,q,t) because D(a, q, oo) = 0.) The real 
exponential function in Eq. ( 1T91) for F(q,t) appears also in Eq. (T4"5l) for D(a,q,t), and this 
function is well behaved at all T. Eqs. ( TT9l) and ( H5l) for the time correlation functions are 
computable as triple iV-fold sums. 
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